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Abstract 

A quantum Monte-Carlo is proposed to describe fusion/fission processes 
when fluctuation and dissipation, with memory effects, are important. The 
new theory is illustrated for systems with inverted harmonic potentials coupled 
to a heat-bath. 

To understand fusion and fission processes, few relevant degrees of freedom (DoF) 
that contain important aspects of the nuclear dynamics: relative distance, mass 
asymmetry, deformation are selected. To properly treat these mechanisms, rel- 
evant DoF should be seen as an Open Quantum System (OQS) coupled to the 
surrounding complex environment of irrelevant DoF (called hereafter environment). 
Most recent advances in OQS should be used in these cases, since close to the barrier 
dissipation, fluctuation and memory (also called Non-Markovian) effects are much 
involved in the dynamics. Recently, we have proposed a new method to describe a 
system (S) + environment (E) with a total Hamiltonian given by 

H = hs + hE + Q®B, (1) 

where hs and He denote the system and environment Hamiltonians respectively 
while Q ® -B is responsible for the coupling. The system density evolution, denoted 
by ps, can be exactly replaced by a stochastic integro-differential master equation 
written as: 

dt /■* 

dps = -J- [hs, Ps] + dt[Q, Ps] dsD{t - s) {Qis))^ 

+ dat)[Q,Ps] + dvit){Q-{Q)s,Ps} (2) 

with 

d^{t) = dt Di{t — s)dvE{s) — dt / D{t — s)duE{s) — idvs{t), drj{t) = dus{t). 
Jo Jo 

Noises dus/E and dvs/E are Gaussian stochastic variables with zero mean values and 
variances equal to: 



dusduE = dvsdvE = ||, dusdvE = dvgduE = 0. (3) 

In Eq. (I2D, D and Di denote time memory functions characterizing the environment 
[2]. 
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To illustrate the method, we consider the Caldeira-Leggett (CL) model which can be 
analytically solved and consists of a harmonic system coupled to a bath of oscillators 
[3] with: 

= + 1? + = ? (it: + r^-"^ 

while B = —J2n'^nXn [!]• Here, he = J2n ^n/ i'^^n^n) is the counter-term that 
insures that the physical frequency is ujo- The thermal behavior of the environment 

is completely defined by its spectral density [1,5] J(uj) = — 6(uj — Un)- In 

the following, a Drude spectral density [2] 

2Mr/ 

= — — ^ ,2 I 02 ' 

is considered where M is the nucleon mass, and rj denotes the strength of the cou- 
pling. In this model, the master equation can be solved by following the stochas- 
tic evolution of first and second moments of P and Q. The QMC efficiency has been 
systematically investigated for various temperatures and coupling strengths in the 
harmonic case [e = 1) [7]. In all cases, averaged evolutions could almost not be dis- 
tinguished from the exact evolution. This is illustrated in figure [T] where the average 
evolution of momentum dispersion, denoted by Spp (left, filled circle), is displayed 
as a function of time and compared to the exact solution (solid line). This approach 
is as accurate as the fourth order in coupling (filled squares) of the Time Convo- 
lutionLess (TCL) method (right). This expansion method consists in a time local 
projective method [1] and leads to master equation with time- dependent transport 
coefficients: 

^J^Psit)= - t[Hs.ps{t)]-'-/\{t)\Q,{Q,ps{t)}]-2t\{t)[Q,{P,Psm 

- ^^[QAQ.Psm+^^^iQAP.psm- (6) 

This method, when second order in the coupling constant is retained, is the preferred 
tool to treat non-Markovian effects in nuclear physics. 

The inverted harmonic oscillator {e = —1) is the first step towards realistic applica- 
tions to fusion and fission. In these physical situations, the passing probability, i.e. 
the probability to pass through the barrier, is of particular importance. It is usually 
defined as: 

P(+oo) = lim-erfc( J^SL ) , (7) 



where q{t) and Sgg(t) denote the expectation value and second moment of Q de- 
duced from the considered theory. We have systematically investigated the difference 
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Figure 1: Evolution of Spp obtained with quantum Monte-Carlo (QMC, left) and 
TCL (right) formalisms are shown as a function of time and compared to the exact 
evolution (solid line). In the TCL case, the second order (TCL2, open squares) or 
fourth order (TCL4, filled squares) in the coupling constant are displayed. Param- 
eters of the simulation are kbT = huQ, ?7 = 0, 5hujQ, hQ = 5hujQ and hujQ = lAMeV. 

between the estimated passing probability and the exact one using the parameter 
AP/P: 

AP ^ P(+00) - Pex( + 00) , . 

P - Pex( + 00) ^ ' 

where P{+oo) and Pea;(+oo) denote the results of the specific calculation considered 
and the exact one respectively. 

Figure [2] presents the evolution of AP/ P as a function of AP, the total system energy 
minus the barrier hight, obtained for different temperatures for the quantum Monte- 
Carlo (filled circles), second order of the TCL method (TCL2, open squares), TCL4 
(filled squares) cases. In this figure, the Markovian approximation is also shown by 
open crosses. 

The TCL4 and the quantum Monte-Carlo methods are in perfect agreement with 
the exact solution for any input parameters. Small differences sometimes observed 
between the stochastic approach and the exact value come from the limited number 
of trajectories used to obtain figure [21 At low AP, TCL2 turns out to be a rather 
poor approximation which confirms the difficulty of treating non-Markovian effects 
below the barrier. It is worth finally to mention that below the barrier, the Marko- 
vian limit gives an even better result than the TCL2 case. 

Quantum Monte-Carlo method is a competitive theory to simulate non-Markovian 
effects in Open Quantum Systems. The present application to inverted parabola is 
the first step towards realistic treatment of fusion/fission. 
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Figure 2: Evolution of AP/P as a function of AE calculated using quantum Monte- 
Carlo (filled circles), TCL2 (open squares), TCL4 (filled squares) and Markovian 
approximation (open crosses) for different temperatures: from left to right T = 5huJ, 
huj and 0, Ihuj (adapted from [7]). 
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